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Molecular phylogenetic relationship of African swallowtail butterflies, Papilio dardanus, P. phorcas, 
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Abstract Molecular phylogenetic relationships of African swallow-tail butterflies Papilio dardanus, P. phorcas, P. 
nireus, and P. demodocus to Japanese papilionid butterflies were analyzed through DNA sequences of mitochondrial 
NDS, COI, and COI genes and are discussed. In addition to 10 butterfly species sequence data obtained in this study, 
other sequence data of the ND5, COI, and COI genes from 20 Papilio species was taken from database (NCBI) to 
construct a phylogenetic tree of Papilio. In particular, we have paid attention to the origin of Batesian mimicry of 
butterflies in both African and oriental (including Japan) regions. The main points of our results are as follows. (1) 
The African mimetic butterfly species P. dardanus and P. phorcas showed a sister relationship, but have no closely 
related species of the genus Papilio in Japan. (2) There was no Japanese species with a clear relationship to the 
African butterfly P. nireus, even though it showed a sister relationship to P. oribazus and P. epiphorbas, which are 
found in Madagascar. (3) The African butterfly P. demodocus showed a close relationship to an oriental species, P 
demoleus with strong support, suggesting that these two species recently diverged from a common origin. This result 
is consistent with that analyzed by use of the COI, COII, and EF-1a genes (Zakharov et al., 2004). (4) The 
phylogenetic relationships of the mimetic species P. protenor, P. macilentus, P. polytes, and P. memnon were resolved 
well, a different result from that reached with lower support by Zakharov ef al. (2004). This result suggests that the 
mimicry in both sexes of P. protenor and P. macilentus is closely related to the female-limited mimicry of P. polytes 
and P memnon. (5) A sister relationship between P. glaucus and P. clytia was suggested with strong support. 

Key words African and Japanese papilionid butterflies, female-limited Batesian mimicry, mimetic butterflies, 
mitochondrial DNA, molecular phylogeny. 


Introduction 


The phylogenetic relationships of swallowtail butterflies of the 
genus Papilio have been widely studied to solve biological 
problems in evolutionary biology, ecology, and genetics. For 
example, the phylogenetic relationships of single species groups 
of three African swallowtail butterflies, Papilio dardanus, P. 
phorcas and P. constantinus were studied to find the origin and 
evolution of the mimetic patterns of P. dardanus. Analysis 
methods have ranged from a cladistic analysis (e.g., Vane- 
Wright and Smith, 1991), hybridization by hand-pairing (e.g., 
Ae, 1979; Clarke ef al., 1991; Sekimura and Ae, 2003) to DNA 
sequencing (e.g., Vane-Wright ef al., 1999 ; Caterino and 
Sperling, 1999; Aubert et al., 1999; Morinaka et al., 1999; Yagi 
et al., 1999, 2006 ; Omoto et al., 2004; Zakharov et al., 2004, 
2008; Simonsen et al., 2011). 


Our main research interest in this paper is to map molecular 
phylogenetic relationships of the African swallowtail butterflies, 


P. dardanus, P. phorcas, P. nireus, and P. demodocus (see Fig. 1) 
to oriental and Japanese papilionid butterflies such as P. xuthus, 
P. polytes, and so on. In particular, we are interested in the origin 


and relationship between mimetic (i.e., Batesian mimicry) 
butterflies in both African and oriental regions including Japan. 
P. dardanus is a species of female-limited mimetic swallowtail 
butterfly widely distributed across sub-Saharan Africa. P. 
dardanus is well known for its spectacular phenotypic 


polymorphism in females (Clarke and Sheppard, 1959 , 1960). 
The females of different geographic races have evolved more 
than a dozen different wing color patterns and have come to 
mimic different species of unpalatable butterflies and moths in 
their specific regions. The males, on the other hand, are 
monomorphic and strikingly different from the females, 
exhibiting a characteristic yellow and black color pattern and 
tailed hind wing. P. phorcas, called the apple-green swallowtail, 
is found in Africa and has been shown to be a sister species of P. 
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Papilio dardanus (male) 





P. nireus 





P. demodocus 
Fig. 1. African swallowtail butterflies, (top left (Papilio dardanus (male), (top right) Papilio 


phorcas (male), (bottom left) Papilio nireus, (bottom right) Papilio demodocus. 


dardanus through some different methods of analysis (e.g., 
Vane-Wright and Smith, 1991; Clarke et al., 1991; Vane-Wright 
et al., 1999). P. nireus, called the narrow green banded 
swallowtail, is also distributed across most of sub-Saharan 
Africa, but is absent in Madagascar and the Comoros. P. 
demodocus, called the citrus butterfly, is distributed across most 
of sub-Saharan Africa including Madagascar and south-west 
Arabia. P. demodocus does remarkably resemble an oriental 
species P. demoleus, which is widely distributed from north-east 
Arabia to India, south-east Asia including Indonesia, the 
Philippines, and also Australia. P. xuthus, called the Asian 
Swallowtail, is one of the commonest swallowtail butterflies in 
Japan, and is widely distributed in central and eastern Asia 
including Russia, China, North- and South Korea. P polytes is a 
female-limited mimetic swallowtail butterfly species widely 
distributed across India and Southeast Asia, including the 
southeast of China, the Philippines, Taiwan, and the Ryukyu 
Islands of Japan (Clarke and Sheppard, 1972). The females of P 
polytes are polymorphic, whereas the male is monomorphic and 
exhibits a white bar on the black hindwing. 


Molecular phylogenetic relationships within the genus Papilio 
have in the past been analyzed by using the DNA sequences of a 
variety of genes (e.g., cyt b, EF-1a, ND5, COI and COII). In 


particular, Yagi ef al. (1999) analyzed the molecular 
phylogenetic relationships of Japanese papilionid butterflies by 
using DNA sequences of the mitochondrial ND5 gene. In the 
current paper, we analyze the phylogenetic relationships of the 
African swallowtail butterflies, P. dardanus, P. phorcas, P. 
nireus, and P. demodocus to oriental and Japanese papilionid 
butterflies by using DNA sequences of mitochondrial ND5, COI, 
COII genes. In the last section, we summarize obtained results 
and examine them in comparison with other results (e.g., Yagi ef 
al., 1999, 2006; Zakharov et al., 2004; Simonsen et al., 2011), 
and we finally discuss the origin of female-limited mimicries of 
both African and Japanese butterflies from other viewpoints than 
molecular phylogeny. 


Materials and methods 
Samples 


The following 10 butterfly species were sequenced for NDS and/ 
or COI-COII genes in this study. 


Papilio protenor (Kasugai, Aichi Pref., Japan, ex T. Sekimura) 


Papilio polytes (Ishigaki-jima Is., Okinawa Pref., Japan, ex Yuji 
Irino) 
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Papilio machaon (Kasugai, Aichi Pref., Japan, ex T. Sekimura) 
Papilio xuthus (Kasugai, Aichi Pref., Japan, ex T. Sekimura) 


Papilio dardanus (Nairobi, Republic of Kenya, ex Steve 
Collins) 


Papilio demodocus (Nairobi, Republic of Kenya, ex Prof. 
Shigeru Ae via Steve Collins) 


Papilio nireus (Nairobi, Republic of Kenya, ex Prof. Shigeru Ae 
via Steve Collins) 


Papilio phorcas (Nairobi, Republic of Kenya, ex Prof. Shigeru 
Ae via Steve Collins) 


Graphium sarpedon (Kasugai, Aichi Pref., Japan, ex T. 
Sekimura) 


Pachliopta aristolochiae (Ishigaki-jima Is., Okinawa Pref., 
Japan, ex Yuji Irino) 


The nucleotide sequence data reported in this study have been 
deposited in the DDBJ, EMBL and GenBank nucleotide 
sequence databases under the accession numbers (NDS: 
LC 014936 -LC 014943 and LC 189140 , COI: LC 189126 - 
LC 189132 , COII: LC 189133 -LC 189139). In addition to the 
sequence data obtained from this study, other sequence data of 
ND5, COI, and COII genes from 20 Papilio specimens were 
taken from GenBank database (NCBI) to construct a 
phylogenetic tree of Papilio (Caterino and Sperling, 1999 ; 
Morinaka et al., 1999 ; Omoto et al., 2004 ; Simonsen et al., 
2011 ; Yagi et al., 1999 ; Yagi et al., 2006 ; Zakharov and 
Hellmann, 2008; Zakharov et al., 2004). 


DNA extraction, PCR and sequencing 


Total DNA was prepared from thorax muscle using a QIAamp 
DNA Mini Kit (QIAGEN, Germany). DNA for each specimen 
was finally dissolved in 200 ul elution buffer, and 2 wl of the 
DNA solution was used as a template for PCR amplification of 
mitochondrial DNA fragments of ND5, COT, and COI genes. 
The PCR primers were according to Su ef al. (1996) for ND5 
gene, Su et al. (2004) for COI gene, and Caterino and Sperling 
(1999) for COI gene. Direct sequencing of PCR products was 
performed with an automated ABI PRISM 3100 Genetic 
Analyzer. 


Phylogenetic Analysis 


The nucleotide gene sequences were aligned by MAFFT (Katoh, 
et al., 2005) with default settings. For phylogenetic analysis, 
model testing was conducted for the sequence alignment with 
MEGA6 (Tamura et al., 2013) under the BIC and AICc criteria, 
and both the criteria selected the GTR+G+I model as the best-fit 
model. Phylogenetic analysis was also performed on MEGA6 
(Tamura et al., 2013) with maximum likelihood (ML) method 


under the GTR+G+I model and 1000 bootstrap replications. The 
phylogenetic tree was outgroup-rooted using the two non-Papilio 
species (Graphium sarpedon, Pachliopta aristolochiae) . 


Results and Discussion 
Phylogenetic relationships of Papilio species 


Mitochondrial ND5, COI, and COII gene sequences were used 
to investigate the phylogenetic relationships of Papilio species, 
especially the position of four African species in relation to 
oriental (in particular, including Japan) species. To compare 
with the results of previous studies, the representative Papilio 
species of the lineages found in Zakharov et al. (2004) were 
included into our phylogenetic analyses. The final sequence 
alignment was 2180 nucleotides in length (ND5: 748; COI: 709; 
COI. 723). 


All the nucleotides of codons (1ST+2nd+3rd) were used for 
phylogenetic analysis, and the resulting tree is shown in Fig. 2. 
To investigate the phylogenetic position of an African Papilio 
species, P. nireus, which is not analyzed in Zakharov et al. 
(2004), we also constructed a phylogenetic tree (Fig. 3) using 
only COI-COII gene sequences, because we did not have NDS 
sequence data of its related species (P. epiphorbas and P. 
oribazus). When compared to the results from Zakharov ef al. 
(2004), there were three main findings obtained from our 
analyses. First, the phylogenetic relationships of five Papilio 
species (P protenor, P. macilentus, P. polytes, P nemnon, and P. 
helenus) clustered in clade C were well resolved with strong or 
moderate support (Figs 2 and 3). The topology was ((((P. 
protenor, P. macilentus), P. polytes) , P. memnon), P. helenus), 
which is different from the result with low support by Zakharov 
et al. (2004) . Our result suggests that the both-sexes mimicry of 
P. protenor and P. macilentus (e.g., Zakharov et al., 2004), 
which are widely assumed to mimic an unpalatable butterfly 
species Byasa alcinous, is closely related to the female-limited 
mimicry of P. polytes and P. memnon (e.g., Clarke and Sheppard, 
1972; Clarke et al., 1968) . Second, P. nireus was clustered with 
P. epiphorbas and P. oribazus with robust support (clade H in 
Fig. 3), supporting the morphological view that these three 
species belong in the same species group (see Zakharov et al., 
2004) . However, the relationships among the three species were 
unclear due to the low support (Fig. 3). Third, the sister 
relationship between P. glaucus (clade J) and P clytia (clade K) 
was supported with high bootstrap value (Figs 2 and 3), 
suggesting that the two species may belong to a common clade. 


Phylogenetic relationships of four African swallowtail butterflies, 
Papilio dardanus, P. phorcas, P. nireus, and P. demodocus to 
Japanese papilionid butterflies were analyzed in the current 
study. The African mimetic species P. dardanus was clustered 
with P. phorcas with strong support (clade A in Fig. 2), but has 
no closely related species of the genus Papilio in Japan (Figs 2 
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Fig. 2. Phylogenetic tree based on mitochondrial ND5, COI and COI gene sequences of 


Papilio species. The tree was constructed by the Maximum Likelihood method based 


on GTR+G+I model. Bootstrap values are shown at nodes, and those less than 50% are 


omitted. The scientific names of species analyzed in this study are shown in boldface. The 


symbol of lineages is according to Zakharov et al. (2004). 


and 3). There was no Japanese species that is clearly closely 
related to the African butterfly P. nireus, even though it showed a 
sister relationship to P. oribazus and P. epiphorbas, which are 
found in Madagascar. The African butterfly P. demodocus 
showed a close relationship to an oriental species, P. demoleus 
with strong support, suggesting that these two species recently 
diverged from a common ancestor. These results are consistent 
with those analyzed by use of the COZ, COII, and EF-la genes 
(Zakharov et al., 2004). 


Female-limited mimicries in the African and oriental regions 


It is known that some species of Papilio exhibit mimetic wing 
patterns world-wide (e.g., Zakharov ef al., 2004 ; d'Abrera, 
1981), and one of their most striking features is female-limited 
mimetic polymorphism. Females of the African P. dardanus 
have more than a dozen different wing patterns or forms, which 
mimic wing patterns of model butterflies and moths, while the 
male has only a simple black and white colored wing pattern 
(Clarke and Sheppard, 1959, 1960). We showed that P. 
dardanus has no closely related species of the genus Papilio 
including P. polytes and P. memnon, which are found in Japan. 
Both the oriental species P. polytes and P. memnon are also 


known to show female-limited mimetic polymorphism even 
though numbers of different female forms are much smaller than 
seen in P. dardanus. Females of P. polytes mimic unpalatable 
butterflies such as Pachliopta aristolochiae and Atrophaneura 
polyeuctes, while females having tails on the hind wings of P. 
memnon are assumed to mimic the unpalatable Atrophaneura 
coon. A close relationship of P. polytes and P. memnon was 
suggested previously by analysis of the ND5 gene (Yagi et al., 
1999). On the other hand, Zakharov et al. (2004) suggested that 
P. polytes is closely related to P. protenor, but not to P. memnon 
on the basis of different genes (the COI, COII, and EF-1a 
genes). The current analysis by a combination of three genes 
(the ND5, COI, and COI) has shown that P memnon seems to 
be a sister taxon to P polytes-P. macilentus-P. protenor, (Fig. 2), 
a result which is different from those of both Yagi et al. (1999) 
and Zakharov et al. (2004). We note here that our current 
analyses and results are based on information from DNA 
sequences of mitochondrial genes. A definitive and correct 
solution to this kind of subtle problem will require further more 
extensive analyses, which will be the task of future research. The 
above-mentioned two different kinds of features (i.e., mimicry 
and close phylogenetic relationship) have been also reported in 
the phylogenetic analysis by Zakharov et al. (2004). The 
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Fig. 3. Phylogenetic tree based on mitochondrial CO/ and COI gene 


sequences of Papilio species. For details, see Fig. 2. 


coincidence between both research results is not accidental, but 
proves its validity, and it seems to show a cue to reveal the 
mechanism of female-limited mimetic polymorphism. Our result 
in Fig. 2 also indicates that the both-sexes mimetic butterflies P 
protenor and P. macilentus are closely related to female-limited 
mimetic butterflies P. polytes and P. memnon and this might 
suggest a hint to approach to the mechanism of both-sexes 
mimicry of P. protenor and P. macilentus. 


In addition to the phylogenetic study, direct research for genes 
controlling mimicry has recently been done intensively and 
enthusiastically. Using bacterial artificial chromosome (BAC) 
sequencing, recombination analyses and comparative genomics, 
Timmermans et al. (2014) showed that the polyallelic H locus 
controlling a set of distinct mimetic phenotypes of P. dardanus 
co-segregates with an interval of less than 500 kb that 1s 
collinear with two other Lepidoptera genomes and contains 24 
genes, including the transcription factor genes engrailed (en) 
and invected (inv). Using an integrative approach combining 
genetic and association mapping, transcriptome and genome 
sequencing, and gene expression analyses, Kunte et al. (2014) 
showed that a single gene, doublesex (dsx) , controls supergene 


mimicry in P. polytes. Further, Nishikawa et al. (2015) showed 
the detailed molecular mechanism of the mimicry in P. polytes, 
where a single 130-kb autosomal inversion, including dsx, 
between mimetic (H-type) and non-mimetic (h-type) 
chromosomes in P. polytes is associated with the mimicry-related 
locus H. They finally proposed the molecular mechanism that 
dsx(H) switches the coloration of predetermined wing patterns 
and that female-limited polymorphism is tightly maintained by 
chromosomal inversion. 
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摘 X 


ミト コン ドリ ア 遺 伝 子 に よる 4 種 の アフ リカ 産 ア ゲハ チョ 
ウゥ の 系 統 関係 AAR ADDA + BREED) 


アフ リカ 産 の アゲ ハチ ョ ウ : オス ジロ アゲ ハ (Papilio 
dardanus), フォ ルカ スミ ドリ アゲ ハ (Papilio phorcas), 

ルリ アゲ ハ (Papilio nireus), アフ リカ オナ シア ゲハ (Papilio 
demodocus) と 日 本 に 生息 する アゲ ハチ ョ ウ 科 の チョ ウ と 
の 分 子 系 統 関 係 を . ミト コン ドリ ア DNA の NDS, Core 
COM 遺伝 子 領域 の 塩基 配列 を 使っ て 分 析 し , 検討 し た . 

アゲ ハチ ョ ウ 属 (Papilio) の 分 子 系 統 樹 を 作成 する に あ 
た っ て は , 本 研究 で 得 ら れ た 10 種 の チョ ウ の 塩基 配列 デー 
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タ の 他 に 20 種類 の ND5, COI と CO 遺伝 子 デ ー タ を 米 
国 国立 生物 工学 情報 セン ター (NCBI) か ら 取 得 し て 利用 
し た . 特に , 本 研究 で は アフ リカ と 東 ・ 南 アジ ア ) 日 本 を 
含む ) 地域 に 生息 する チョ ウ の ベー ツ 擬 態 の 起源 に 注目 し 
た . 本 研究 で 得 ら れ た 主 な 結果 は 以下 の 通り で ある . (1) 

オス ジロ アゲ ハ (P dardanus) と フォ ルカ スミ ドリ アゲ ハ 
(P phorcas) は , 姉妹 関係 を が し た が , 日 本 に 生息 する ア 
ゲハ チョ ウ 属 の チョ ウ と 分 子 系 統 的 類縁 関 作 を 持た な い . 

(2) ルリ アゲ ハ (P nireus) lk, マダ ガス カル に 生息 する 
オリ バズ スル リア ゲハ (Poribazus) そし て エピ フォ ル バ 
スル リア ゲハ (P epiphorbas) と 姉妹 関係 を 示し た . し か し , 
日 本 産 の アゲ ハチ ョ ウ 類 と は 上 類縁 関係 を 示さ な か っ た . (3) 
アフ リカ オナ シア ゲハ (P. demodocus) tk, アジ ア の 東部 お 
よび 南部 地域 に 広く 生息 する オナ シア ゲハ (の demoleus) 
と 姉妹 関係 を 示 した. これ は , この 2 種 の チョ ウ が 近 緑 で 











あり 共通 祖先 由来 で ある こと を 示唆 し て いる . この 結果 は 
COI, COII と EF-Ia 遺伝 子 を 用 いた Zakharov et al. (2004) 
の も の と 一 致す る . (4) クロ アゲ ハ (P protenor) と オナ 
ガ ア ゲハ (P. macilentus) に 観察 され る 典 雄 両性 擬態 は シ 
ロ オ ビア ゲハ (P polyte) と ナガ サキ アゲ ハ (P memnon) 
に 見 られ る 雌 限定 の 擬態 と 緊 窒 な 関係 が ある こと が 示唆 さ 
れ た . (5) P glaucus (clade J) & P clytia (clade K) の 姉 
妹 関 係 は 強く 支持 され た . 最後 に . アフ リカ と 東 ・ 南 アジ 
Y (日 本 を 含む ) に 生息 する 蝶 に 見 られ る 雌 限 定 擬態 に つ 
いて , 分 子 系 統 学 的 議論 だ け で な く 他 の 視点 (例え ば , 最 
近 発見 され た 目的 遺伝 子 ) か ら も 議論 した. 
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